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Qu Abstract 

In this paper we present the C++ package CRunDec which implements all relevant formulae 
needed for the running and decoupling for the strong coupling constant and light quark 
masses. Furthermore, several formulae are implemented which can be used to transform the 
heavy quark masses among different renormalization schemes. CRunDec is the C++ version 
on the Mathematica package RunDec containing several updates and improvements. 

PACS numbers: 12.38.-t 12.38.Bx 14.65.-q 
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Program summary 

Title of program: CRunDec 
Available from: 

[http : //www-ttp . physik . uni-karlsruhe . de/Progdata/ttpl2/ttpl2-02/ 

Computer for which the program is designed and others on which it is operable: Any 
computer where a C++ compiler is running. 

Operating system or monitor under which the program has been tested: Linux, Windows 
No. of bytes in distributed program including test data etc.: 85 000 
Distribution format: source code 

Keywords: Quantum Chromodynamics, running coupling constant, running quark 
mass, on-shell mass, MS mass, decoupling of heavy particles 

Nature of physical problem: The value for the coupling constant of Quantum Chro- 
modynamics, ai n/ ^(/i), depends on the considered energy scale, /i, and the number of 
active quark flavours, n/. The same applies to light quark masses, m^^/i), if they are, 
e.g., evaluated in the MS scheme. In the program CRunDec all relevant formulae are 
collected and various procedures are provided which allow for a convenient evaluation 
of a,s\fJ>) and rriq lf \fj,) using the state-of-the-art correction terms. 

Method of solution: CRunDec is implemented in C++. For the solution of the differential 
equations an adaptive Runge-Kutta procedure has been implemented. 

Restrictions on the complexity of the problem: It could be that for an unphysical choice 
of the input parameters the results are nonsensical. 

Typical running time: In general the run time for the individual operations is below a 
millisecond. 



1. Introduction 

Among the fundamental quantities of Quantum Chromodynamics there are certainly 
the anomalous dimensions which control the scale dependence of the parameters and fields. 
In this context a particular role is taken over by the beta function and anomalous mass 
dimension in the MS scheme which govern the running of the strong coupling constant a s 
and the quark masses m q . Every time a flavour threshold is crossed in the running process 
decoupling relations have to be applied which guarantee that the heavy quark is integrated 
out from the theory. More than ten years ago the Mathematica package RunDec has been 
written [lj which incorporates all relevant formulae. However, for many application it is 
more convenient to have the running and decoupling routines available in the framework of 
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a commonly used programming language and not within a computer algebra system. Thus 
we have decided to convert the most important routines to C++ adding at the same time new 
routines and improving others w.r.t. the Mathematica version. In what follows we describe 
the C++ file CRunDec which constitutes a C++ class and contains the functions known from 
RunDec as public components. 

The use of CRunDec does not require any knowledge about object oriented programming. 
The following skeletons exemplify the usage and can easily be adapted to the problem at 
hand. It is possible to work with pointers to an object of the type CRunDec and access the 
member functions correspondingly: 

#include <iostream> 
#include "CRunDec. h" 
using namespace std; 
int main(){ 

CRunDec * <pointer> = new CRunDecO; 

double <result> = <pointer> -> <function>(<parameters>) ; 
return (0) ; 

} 

Alternatively also the following realization is possible: 

#include <iostream> 
#include "CRunDec. h" 
using namespace std; 
int main(){ 

CRunDec <object>; 

double <result> = <object> . <function>(<parameters>) ; 
return (0) ; 

> 

Explicit examples will be given below. 

The remainder of the paper is organized as follows: In the next Section all available 
functions are described, Section [3] contains typical examples which exemplify the usage of 
CRunDec. 

2. Structure of CRunDec 

All public components of the C++ class CRunDec consist of functions which have the same 
name as the corresponding function in the Mathematica version In the following we list 
the function header (including the variable names; see also Tab. [1]) which — together with 
the description in the Appendix of Ref. [l[ - specifies both the usage and the purpose of 
the function. There are some additions implemented in CRunDec which are also described. 

Let us mention that all functions listed in Subsections I2.2ff2.4l are overloaded w.r.t. the 
argument rif (number of active flavours). I.e., this argument can be omitted in case n/ is 
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symbol m C++ code 


mathematical symbol 


meaning 


asmu, asmuO, asmul 


OisW, ot 8 {fi Q ), 


strong coupling constant 


mq , mqO 




MS quark mass 


mUo , mMS 


M q , m q 


on-shell and MS quark mass 


mRI 


m RI 


regularization invariant mass 


mRGI 


m 


renormalization group invariant mass 


Mth 


Ma, 


heavy quark mass 


mu, muO, mul 


A*, Mo, A*i 


renormalization scale 


muth 


Mth 


decoupling scale 


nf 




number of active flavours 


nloops 


number of loops 



Table 1: Meaning of the variables used in the function headers. 



specified in the declaration of the CRunDec object (see also examples in Section [3]). In this 
context two auxiliary functions are quite useful: GetNf () returns the specified number of 
active flavours and SetNf (int nf ) can be used to set the number of active flavours. 

2. 1 . Input parameters 

For convenience of the user some frequently used input parameters arepre-defined in the 
file CRunDec .h and can be used during the calculation. They read [1, 0, 0, 0] 

#define asMz 0.1183 
#define Mz 91.18 
#define Mt 173.2 
#define Mb 4.8 
#define Mc 1.5 
#define muc 1.279 
#define mub 4.163 
#define Mtau 1.777 

If not stated otherwise these numbers are used in the examples presented in Section [3j In 
case other numerical values shall be used it is straightforward to redefine the preprocessor 
variables. The examples given in the appendix of Ref. [l| can be reproduced with 

#define asMz 0.118 
#define Mz 91.18 
#define Mt 175. 
#define Mb 4.7 
#define Mc 1.6 
#def ine muc 1 . 2 
#define mub 3.97 
#define Mtau 1.777 
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2.2. Functions related to the running of a s and m q 

• double LamExpl (double asmu, double mu, int nf, int nloops) ; 

• double Lamlmpl (double asmu, double mu, int nf, int nloops); 

• double AlphasLam (double Lambda, double mu, int nf, int nloops); 

• double AlphasExact (double asmuO, double muO, double mul, 

int nf, int nloops); 

• double mMS2mMS (double mqO, double asmuO, double asmul, 

int nf, int nloops); 

• AsmMS AsmMSrunexact (double mqO, double asmuO, double muO, 

double mu, int nf, int nloops); 

The function AsmMSrunexact is new in CRunDec. It solves simultaneously the differential 
equations for a s and m q with initial values m q (no) and a s (/xo) and rif active quark flavours. 
The return type of AsmMSrunexact is a struct with two double components, a s (fi) and m q (ii). 
The corresponding code in CRunDec looks as follows: 

struct AsmMS { 

double Asexact; 
double mMSexact ; 

>; 

For convenience of the user there is a pre-defined variable AsmMS AM. Both components of 
AM are initialized to zero when creating a CRunDec object. 

2.3. Functions relating different mass definitions 

• double mOS2mMS (double mOS, double mq[] , double asmu, double mu, 

int nf, int nloops); 

• double mMS2mOS (double mMS, double mq[] , double asmu, double mu, 

int nf, int nloops); 

• double mOS2mMSrun (double mOS, double mq[] , double asmu, double mu, 

int nf,int nloops); 

• double mMS2mOSrun (double mMS, double mq[] , double asmu, double mu, 

int nf,int nloops); 

• double mOS2mMSit (double mOS, double mq[] , double asmu, double mu, 

int nf,int nloops); 

• double mOS2mSI (double mOS, double mq[] , double asM, 

int nf, int nloops); 
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• double mMS2mSI (double mMS, double asmu, double mu, 
int nf, int nloops) ; 



• double mMS2mRI (double mMS, double asmu, int nf, int nloops); 

• double mRI2mMS (double mRI, double asmu, int nf , int nloops) ;0 

• double mMS2mRGI (double mMS, double asmu, int nf, int nloops); 

• double mRGI2mMS (double mRGI , double asmu, int nf, int nloops); 

• double mMS2mRGImod (double mMS, double asmu, int nf, int nloops); 

Note that the light-quark-mass effects can be taken into account with the help of the 
array mq [] which is defined as 

double mq[4] ; 

By default all elements of mq[] are zero. In case non-zero values are needed the array has 
to be filled before the corresponding function is called. 

In CRunDec the implementation of mMS2mSI has been modified as compared to the 
Mathematica version. It is now based on AsmMSrunexact and avoids the computation of 
Aqcd in intermediate steps which is perturbatively more stable, in particular for lower renor- 
malization scales. Similar modifications have been performed in mMS2mOSrun. 

The function mMS2mRGImod is new in CRunDec. It is defined in analogy to mMS2mRGI, 
however, the more commonly used convention has been adopted where the function c(x) in 
Eq. (11) of Ref. [l| is evaluated for x = 2(3 a s /7r instead of x = a s /ir. 

2.4- Functions related to the decoupling of heavy quarks 

As compared to the Mathematica version CRunDec contains the decoupling relations only 
for the case of on-shell heavy quarks which are most relevant for the practical purposes. 
Furthermore, the functions DecLambdaUp and DecLambdaDown have not been implemented 
in the C++ version since it is recommended to use A1L2A1H and A1H2A1L in case a flavour 
threshold is crossed during the running of a s . 

• double DecAsDownOS (double asmu, double Mth, double muth, 

int nf, int nloops); 

• double DecAsUpOS (double asmu, double Mth, double muth, 

int nf, int nloops); 

• double DecMqDownOS (double mq, double asmu, double Mth, double muth, 

int nf, int nloops); 



^ote a typo in the example to mRI2mMS in the Appendix of Ref. 1.]: mRI2mMS [175 , . 107 , 175, 6, 3] 
should read mRI2mMS [175 , . 107 , 6 , 3] . 
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• double DecMqUpOS (double mq, double asmu, double Mth, double muth, 

int nf, int nloops) ; 

In CRunDec the functions DecAsDownOS and DecAsUpOS also contain the four-loop decou- 
pling relations which have been computed in Refs. 0, @|. Note that in the functions of this 
subsection the parameters nj refers to the number of flavours in the effective theory [lj]. 

2. 5. Functions related to the combination of running and decoupling 

• double A1L2A1H (double asmuO, double muO, TriplenfMmu decpar[], 

double mul, int nloops); 

• double A1H2A1L (double asmuO, double muO, TriplenfMmu decpar[], 

double mul, int nloops); 

• double mL2mH (double mqO, double asmuO, double muO, 

TriplenfMmu decpar[], double mul, int nloops); 

• double mH2mL (double mqO, double asmuO, double muO, 

TriplenfMmu decpar[], double mul, int nloops); 

The parameters governing the decoupling are contained in the array decpar[] where 
each element contains the triple {rif, Mth, //th} which is realized in the structure 

struct TriplenfMmu { 
int nf ; 
double Mth; 
double muth; 

>; 

There is a pre-defined variable TriplenfMmu nfMmu[4] ; which can be used when calling 
the above functions. Note that the components of decpar are set to zero at the end of the 
above functions. 

In CRunDec we refrain to implement the function AsRunDec which automatically deter- 
mines the number of active flavours for the initial and final energy scale and performs the 
corresponding running and decoupling steps. In practice it turns out that the decoupling of 
the heavy quark with mass Mth at the scale /x t h = M t h is not convenient for all applications. 
Furthermore, there are ambiguities as far as the number of active flavours is concerned in 
case a s (M t h) has to be evaluated using AsRunDec. Thus, it is recommended to use A1L2A1H 
and A1H2A1L instead. 

3. Typical examples 

In this section we present some typical examples which exemplify the usage of CRunDec. 
In the following we only display the part of the code related to CRunDec; the complete 
programs can be found in the file example. cc which comes together with CRunDec. 
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Running of a s from M t to Mb with five active flavours 

It is either possible to create a CRunDec object where five active flavours are already 
specified 

CRunDec* p0bjnf5 = new CRunDec(5); 

or leave rif unspecified 

CRunDec* pObj = new CRunDecO; 

In the former case af\Mb) is computed from af\M t ) = 0.108 via 

p0bjnf5 -> AlphasExact (0 . 108, Mz, Mb ,4) ; 

whereas in the latter case one has 

pObj -> AlphasExact (0. 108, Mz, Mb, 5,4) ; 

In both evaluations four-loop accuracy is assumed leading to the result a^(Mb) = 0.183. 

The number of active flavours for the object pObjnf 5 can be obtained with the help of 
p0bjnf5 -> GetNfO; and pObj -> SetNf(5); sets n f = 5 for pObj . 

Compute af\M z ) from af\M T ) 

Assuming a value of the strong coupling as extracted form r decay as (M T ) = 0.332 [i| 
the task is the computation of af'(Mz)- If the decoupling of the charm and bottom quark 
is performed for /i th = 2M C and /i t h = M b , respectively, one has to specify 

CRunDec crundec; 
crundec .nfMmu[0] .nf = 4; 
crundec .nfMmu[0] .Mth = Mc; 
crundec .nfMmu[0] .muth = 2*Mc; 
crundec .nfMmu[l] .nf = 5; 
crundec .nfMmu[l] .Mth = Mb; 
crundec .nfMmu[l] .muth = Mb; 

Afterwards af'(Mz) is computed via 

crundec . A1L2A1H (0 . 332 , Mtau , crundec . nf Mmu , Mz , 4) ; 

where four-loop accuracy for the running (corresponding to three-loop decoupling relations) 
has been assumed. As a result one obtains af\M z ) = 0.1200. 
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Compute m^nib) from mb(10 GeV) 

In Ref. j^] the MS bottom quark mass has been extracted for /i = 10 GeV as m^lO GeV) = 
3.610 ± 0.016 GeV. The scale-invariant mass mb(mb) can been computed with four-loop ac- 
curacy via 

CRunDec cmpasmb(5); 

double asmu = cmpasmb . AlphasExact (0 . 1189, Mz, 10,4) ; 
cmpasmb . mMS2mSI (3.610, asmu ,10,5,4); 

where af\Mz) = 0.1189 has been used. The result reads m^m^) = 4.163 GeV. In the 
Mathematica version of RunDec it is not recommended to use mMS2mSI since in intermediate 
steps Aqcd is used and thus the final result does not have the required precision. 

Compute m t (mt) from the top quark on-shell mass 

The on-shell top quark mass has been measured at the Tevatron experiments DO and 
CDF to M t = 173.2 ± 0.9 GeV Q. If this value shall be transformed to m t (m t ) one proceeds 
in the following way 

CRunDec cmpmt(6); 
double MtOS = 173.2; 

cmpmt .nf Mmu[0] .nf = 6; 
cmpmt.nfMmu[0] .Mth = MtOS; 
cmpmt .nfMmu[0] .muth = MtOS; 

double as6Mt = cmpmt . AlL2AlH(asMz,Mz, cmpmt .nfMmu, MtOS ,4) ; 
double mtMt = cmpmt .mOS2mMS (MtOS, cmpmt .mq, as6Mt , MtOS ,3) ; 
double mtmt = cmpmt .mMS2mSI (mtMt, as6Mt, MtOS, 4) ; 

The final result reads mt(mt) = 164.0 GeV. 

Compute on-shell charm and bottom quark mass from m c (m c ) and rrib(mb) 

After defining a pointer (p0bjnf5) to a CRunDec object with five active flavours (as in 
the example above) one obtains a^m^mj,)) via 

double alpha5mub = p0bjnf5 -> AlphasExact (asMz, Mz, mub, 4); 

and subsequently the on-shell mass Mb with two and three-loop accuracy with the help of 

double mb0S2 = p0bjnf5 -> mMS2mOS(mub, p0bjnf5->mq, alpha5mub, mub, 2); 
double mb0S3 = p0bjnf5 -> mMS2mOS(mub, p0bjnf5->mq, alpha5mub, mub, 3); 

The results read M b = 4.762 GeV and M b = 4.909 GeV, respectively. 

In the case of the charm quark one proceeds in analogy by evaluating in a first step 

a;( 4) (m c (m c )) 
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pObjnf 5 -> nfMmu[0] .nf =4; 
p0bjnf5 -> nfMmu[0] .Mth = Mc; 
p0bjnf5 -> nfMmu[0] .muth = 2*Mc; 

double alpha4muc = p0bjnf5 -> AlH2AlL(asMz , Mz, p0bjnf5 -> nfMmu, mub, 4); 

and afterwards the on-shell mass M c to two and three loops 

double mcOS2 = p0bjnf5 -> mMS2mOS (muc , pObjnf 5->mq, alpha4muc, muc, 2); 
double mcOS3 = p0bjnf5 -> mMS2mOS(muc, pObjnf 5->mq, alpha4muc, muc, 3); 

leading to M c = 1.494 GeV and M c = 1.573 GeV. 
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